Non-local fluctuation correlations in active gels 
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Many active materials and biological systems are driven far from equilibrium by embedded agents 
that spontaneously generate forces and distort the surrounding material. Probing and characteriz- 
ing these athermal fluctuations is essential for understanding the properties and behaviors of such 
systems. Here we present a mathematical procedure to estimate the local action of force-generating 
agents from the observed fluctuating displacement flelds. The active agents are modeled as oriented 
force dipoles or isotropic compression foci, and the matrix on which they act is assumed to be ei- 
ther a compressible elastic continuum or a coupled network-solvent system. Correlations at a single 
point and between points separated by an arbitrary distance are obtained, giving a total of three 
independent fluctuation modes that can be tested with microrheology experiments. Since oriented 
dipoles and isotropic compression foci give different contributions to these fluctuation modes, ra- 
tiometric analysis allows us characterize the force generators. We also predict and experimentally 
find a high-frequency ballistic regime, arising from individual force generating events in the form 
of the slow build-up of stress followed by rapid but finite decay. Finally, we provide a quantitative 
statistical model to estimate the mean filament tension from these athermal fiuctuations, which 
leads to stiffening of active networks. 

PACS numbers: 87.16.Ka, 62.20.D- 
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I. INTRODUCTION 

Living organisms are inherently non-equilibrium sys- 
tems that continually harvest external resources in order 
to maintain their activity and complex hierarchical or- 
dering P, [^. Thermodynamic equilibrium is therefore 
expected to apply, at best, to length and time regimes 
that are not associated with biological processes. The 
cellular cytoskeleton, a complex assembly composed of 
filamentous proteins, motor proteins and other acces- 
sories which provides a cellular scaffolding [H-Q, is a 
case in point. This dynamic structure hydrolyses en- 
ergy transfer molecules such as adenosine triphosphate 
(ATP) to promote morphological and mechanical varia- 
tions, thus regulating many biomechanical processes such 
as organelle transport, muscle contraction, and cell divi- 
sion [3t3| • Since such processes need not obey the statis- 
tics of thermodynamic equilibrium jSj, |9j , novel analytic 
techniques need to be devised both to predict the vis- 
coelastic response, and to physically describe the active 
entities from experimental data. 

The violation of thermodynamic equilibrium has been 
demonstrated in model in vitro cytoskeletons consisting 
of an actin filament network coupled to the motor pro- 
tein myosin JlOl, which generates ~ pN forces upon ATP 
hydrolysis [3|, The material's mechanical viscoelas- 
tic response was simultaneously determined by two mi- 
crorheology protocols, one (passive) that assumes ther- 
modynamic equilibrium and a second (active) that does 
not [ll| - [l6| . Any difference between these two protocols, 
which has been observed at low [i.e. biological) frequen- 
cies, is attributed to the existence of non-equilibrium 
fiuctuations [l0|, [H, Ezl- A departure from Gaussian 
fluctuation statistics has also been observed in active, en- 



tangled actin-microtubule mixtures |18| . The existence 
of athermal force generators in various model systems 
has been shown to spontaneously create non-equilibrium 
flow patterns and structures [19,-211] , requiring novel ana- 
lytical techniques to understand their collective behavior 
[22h28| . Introducing permanent crosslinkers such as bi- 
otin/avidin complexes into in vitro cytoskeletal systems 
blocks filament flow and reduces complexity to that of 
fluctuations in a static medium. This allows for a physi- 
cal characterization of the athermal fluctuations, leading 
the way to a fundamental understanding of their biolog- 
ical implications. 

Here we derive the spatial and temporal correlations 
in the fluctuating displacement fleld generated by a pop- 
ulation of athermal force generators, or 'firers', dispersed 
in a linear viscoelastic material, in a form that can be di- 
rectly compared to microrheology data. Firers are mod- 
eled as force dipoles, generating stress and strain fields 
that propagate long distances as dictated by material re- 
sponse. As a consequence the correlations induced by 
motor activity are non-local] that is, a single motor acti- 
vation generates correlated motion in distant parts of the 
medium. For the surrounding environment, we consider 
two material types, (i) & single-component elastic con- 
tinuum, and (ii) a, two-component system consisting of a 
sparse elastic network that is frictionally coupled to an in- 
terspersed, incompressible solvent, known as the '2-fluid 
model' |29| - [33 | . The latter is closer in spirit to the 
cellular cytoskeleton, but as our calculations will show, 
the consequence of the explicit solvent-network descrip- 
tion to predicted microrheology data is negligibly small 
in most cases. Thus the single-component model should 
approximately apply to systems of interest, as has been 
already assumed in prior microrheology studies. Since we 
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expect a critical role of strong local compression induced 
by molecular force dipoles on the mechanics of networks, 
we do not assume material incompressibility, in contrast 
to previous approaches @, Q . Macroscopic isotropy and 
homogeneity, relevant to in vitro experiments, are as- 
sumed throughout j33l - l36l |. 

For both classes of material, we derive exact expres- 
sions for the displacement power spectrum Cij{uj), which 
is the magnitude of displacement fluctuations (it^) ex- 
pressed as a function of frequency lo, and the two-point 
spectrum £'ij(R, w) that quantifies correlations in dis- 
placements between points separated by a vector dis- 
tance R. These quantities can be directly compared to 
data from 1 and 2-point microrheology experiments, re- 
spectively. We find 3 sources of variation with frequency: 
(i) The power spectrum of the forces generated by in- 
dividual firers. For an actin-myosin system, in which 
the contractile motor stresses build up gradually but are 
rapidly released, the power spectrum is expected to take 
the form ^ uj^^, as has been observed [13, El, [13, HI] ■ In 
reality the time required to release the developed com- 
pression is short but finite, so to probe this effect we 
solve a simple model with a controllable release time, and 
predict a crossover to a high-frequency ballistic regime 
~ u!^^. Experimental data confirming the existence of 
this ballistic regime is also provided, (ii) The frequency 
dependence of the mechanical response (shear modu- 
lus iJ,{uj)) of the matrix. The linear mechanical response 
of an actin gel has been studied intensively and is now 
relatively well understood as a network of semi-flexible 
polymers [soUi^ . at least for affine deformation (45l - [5^ 
and negligible contribution from crosslinkers or other ac- 
cessory proteins [53 55j. (Hi) When the material is as- 
sumed to obey the 2-fluid model, a third source of varia- 
tion comes from the frictional coupling between the net- 
work and solvent, which becomes tighter at higher fre- 
quencies and suppresses compressional/longitudinal re- 
sponse modes. For parameters relevant to actin-myosin 
systems, our calculations indicate that this contribution 
is expected to be small. 

An example of the correlation functions for the 2-fiuid 
model from our calculations is given in Fig. [l] where we 
have projected out the parallel Z?" (R) and perpendic- 
ular (R) components of Dij (R) in the usual way, as 
Aj(R) = D^R)R^Rj + D^{R){5,j - R,R.j) with the 
Kronecker delta, R = |R| and R = R/i?. These calcu- 
lations are based on the assumption of a population of 
randomly distributed and oriented, independently-firing 
active agents modeled as force dipoles. The inset shows 
the spatial variation of DijiTl) with R (as which is 

consistent with experiments and specific to force dipoles; 
higher-order multipoles give different power laws, as ex- 
plained in Sec. IIII Cl below. The scaling with R does not 
determine the symmetry of the underlying force dipole, 
and we identify two cases of relevance to actin-myosin 
systems, namely an oriented ioice dipole consisting of two 
nearby placed equal-and-opposite forces, corresponding 
to a single myosin minifilament attached to the actin net- 



work, and an isotropic compression center corresponding 
to the orientational average of the first class. This lat- 
ter case is expected to arise due to fluctuations in the 
local concentration of motor protein, or filament rupture 
leading to local network compression. Below we derive 
Cij and Dij for both classes of dipole, explicitly confirm- 
ing Dij ^ 1/R but revealing qualitatively different low- 
frequency behavior for the parallel correlations D^^{R), 
which approaches 2D^{R) for oriented dipoles in the in- 
compressible limit, but entirely vanishes for compression 
centers. Hence, by investigating the length-scale depen- 
dence of each athermal fiuctuation mode or by comparing 
the amplitude of different independent modes, it is now 
possible to identify the class of firer in a system of un- 
known composition, and further test the validity of the 
theoretical modeling or fine-tune the model to specific 
experimental systems. 

Before proceeding, it should be noted that active gels 
are not the only situation where internally-generated 
force impulses are of relevance. Soft glasses and col- 
loidal gels relaxing towards their equilibrium state ex- 
hibit local spontaneous rearrangement events or micro- 
collapses [56|-[58[, generating dipolar force perturbations 
much as in active gels. Elastic propagation of stress from 
local irreversible rearrangements is the essence of various 
physical descriptions of glassy systems (see [59l - l62j and 
references therein), which are also typically modeled as 
a consequence of creation and annihilation of local force 
dipoles. The results obtained in this study thus apply 
to a broader range of non-equilibrium systems than sug- 
gested by the primary, biophysical motivation. 

This paper is arranged as follows. In Sec. [Tl] we de- 
fine both oriented force dipoles and isotropic compres- 
sion centers in terms of a frequency-dependent dipole 
moment k{ui). We show that a simple model with slow 
stress build-up followed by a rapid but finite-duration 
release exhibits a high-frequency ballistic regime in its 
power spectrum, |k(u;)P ~ Microrheology data 

consistent with this prediction is also given. The cor- 
relation functions of the athermal fiuctuating displace- 
ments fields are derived in Sec. IIIIl We first define the 
mechanical response of 2 types of materials, single- and 
two-component, in Sec. IIII Al and then derive Cij{uj) 
and DijCRjio) for oriented force dipoles for both materi- 
als in Sec. UTTB] and Sec. UTTC] respectively. In Sec. HITdI 
the same quantities are derived for isotropic compression 
centers. Finally in Sec. lIVI we estimate the mean network 
tension from the magnitude of athermal fiuctuations in 
an active gel. 



II. INDIVIDUAL FIRERS 

An active agent or 'firer' here refers to any local con- 
struct embedded within a continuum body that is capa- 
ble of spontaneously generating mechanical stress. For 
example, in an actin-myosin system each firer is a bipo- 
lar minifilament of 0(10) myosin motors that consumes 
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FIG. 1: Example of the one-point and two-point correla- 
tion functions Cij{uj) and Dij[R,uj), resp., from the analyt- 
ical expressions (|10p and (|13[l for the two-fluid model. The 
inset shows Dij{uj = ls~^) for different separations R. Fol- 
lowing [ssl ]. the network shear modulus is n{u}) = /io(l -f 
[— ioj/ojo]^''*) with Wo = 10 s~^. For the dipole moment power 
spectrum |K(aj)|^, we used the finite-release model ((3]) with 
T — 100s and a — 10"'^. Other parameters are ly — 0.3 and 
r = lO^'^ kg m~^ s~^ corresponding to a mesh size ^ « lOOnm 
and solvent viscosity i] « 10~^ Pa s. The normalization fac- 
tors are Co = CK^/^i'^a and Do — ckq/iioR with kq — 2t{k)^ . 

the local supply of ATP upon activation [2^ HI, [6^ . 
We consider only the linear response, and leave non- 
linear considerations to future work (such as the mixed 
response mode introduced in [IGj]). The active forces can 
then be written as a multipole expansion of the form 
^ fe"W"5{r), where e, / are microscopic length and 
force scales, resp., and n= 1,2... corresponds to dipoles, 
quadrupoles etc. The monopole term n = vanishes for 
a stable network as the total force must balance, but 
there is no a priori reason to neglect higher terms. Thus 
firers are here taken as dipolar to leading order, and we 
shall later confirm that experiments are in line with this 
assumption. 

A single firer in an active gel at the position x is mod- 
eled as a pair of equal-and-opposite forces ±/n located 
at x=F |n, respectively, where n is the firer orientation, as 
shown in Fig. ^a) . For simplicity we take the limit e ^ 
and / — > cx) with a fixed dipole moment k = ef. This 
reduces each firer to a point force dipole of zero spatial 
extent, which is a valid approximation when addressing 
quantities that vary on length scales much larger than 
the size of the firer (as in the typical microrheology sit- 
uation). The force density due to a single point dipole 
is 

f""'(r, t) = K(t)n(n • V)(5(r - x) (1) 

where the suffix 'ani' differentiates these dipoles, with 
orientation n, from the isotropic compression centers de- 
fined below. The time-dependent strength of the dipole 




(a) (b) 

FIG. 2: (a) A single contractile force dipole of size e, 
force / and orientation n. (b) Schematic representation of 
an isotropic compression center as the orientational average 
of (a). 

moment K{t), which as defined is positive for contractile 
firers and negative for expansive ones, denotes a train of 
activation events, including the possibility of 'dead-time' 
between pulses with «; = 

As mentioned in the introduction, we also consider 
isotropic compression centers as in Fig. [lib), applica- 
ble when local clusters of coupled firers activate near- 
simultaneously or when independent firers are locally ag- 
gregated. Formally these are the average of (jlj over the 
orientation vector fi, 

P°(r,i) = ^«;(0V<5(r-x) (2) 

where we have used the identity hifij — >■ ^Sij un- 
der averaging. Such isotropic force dipoles or dila- 
tion/compression centers have been studied in a number 
of contexts [H, [gI] . 

A. Power spectrum for a single firer 

The temporal dependence of the dipole moment K{t) 
can by itself significantly contribute to the system's fluc- 
tuation spectrum. It is experimentally observed that 
strains induced by myosin minifilaments slowly build up 
over time scales that in some cases extend up to ~1 — 
10 s, and then rapidly release [l3l, allowing the internally 
generated tension to relax to an unstrained state. This 
rapid release generates a pseudo-diffusive power spec- 
trum ^ w"^ for frequencies larger than the inverse char- 
acteristic build-up time, despite the lack of any diffusing 
field [37,_38] ■ At frequencies far below this time, a plateau 
^ w° is expected. However, although rapid, the time for 
stress release is nonetheless finite. As already discussed 
in jlCi], stress build-up is accompanied by significant lo- 
cal network compression deriving from the non-linear re- 
sponse of actin filaments. In order to relax, it is therefore 
necessary for the network to dilate, the rate of which is 
limited by the inflow of solvent. A relaxation time of 
~ 0.1 s was estimated by balancing frictional and elastic 
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FIG. 3: Notation for the finite-release model. The dipole 
moment k increases linearly at a rate kg, initiates release at 
a rate after which it decreases linearly at a rate — kr until 
reaching zero when the cycle restarts. 



forces during relaxation, and found to be consistent with 
experimental observation [l^. To quantify the effect of 
a finite release time on the power spectrum, we consider 
here a simple model in which the dipole moment of a sin- 
gle firer K{t) slowly increases in time before stochastically 
releasing stress as shown in Fig. [3l The model is defined 
by 3 parameters: The rate of increase of the moment k 
in the growth phase, kq > 0; the corresponding rate in 
the release phase, — kr < 0, and the release rate at 
which growth switches to release. Note that k always 
returns to zero before growth restarts, so it is bounded 
and does not diffuse. The steady state solution is found 
following the procedure described in Appendix \^ The 
resulting power spectrum |k(cj)P due to this single firer 
is 
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(3) 



where a = kq/k^ < 1 is the ratio of growth to release 
rates, and (k) = tkq is the mean moment in steady 
state. This expression admits 3 regimes: (i) A low fre- 
quency plateau ~ w° for <C t"-'^; (ii) An in- 
termediate pseudo-diffusive regime |K(a;)p ^ for 
^ a; ^ (ar)"^, where the upper crossover frequency 
is the inverse mean release time; and (in) A high fre- 



quency ballistic regime |K(a;)|^ 



for (jj ^ (ar) 



In the rapid release limit a — >■ 0, the ballistic regime van- 
ishes and the pseudo-diffusive regime extends to all high 
frequencies, as expected. An example for a — 1/100 is 
given in Fig. 21 

For an actin-myosin system, the ballistic regime should 
be clearly resolved for frequencies greater than ~10 Hz 
using the estimated release time ^ 0.1 s discussed above. 
This is around the freq uency where thermal fluctuations 
start to dominate [3 Ell, so achieving a satisfactory 
signal-to-noise ratio is already challenging. Further- 
more this frequency is above the elastic plateau for typ- 
ical network densities, so the ballistic spectrum may be 
obscured by additional frequency variation coming from 
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FIG. 4: The thick line gives the power spectrum of a 

single dipole according to @, with a — kc/kn ~ 1/100, so 
the mean release time is 1% of the growth time. The straight 
dashed hnes have given slopes and both axes have been made 
dimensionless by suitable normahzation. 



the network viscoelasticity. For highly crosslinked net- 
works, however, thermal noise is suppressed and the elas- 
tic plateau extends to higher frequencies [i^ , so the con- 
tribution of the fluctuation of firers, |K(a;)p, should dom- 
inate the frequency dependence of the experimentally- 
measured displacement power spectrum, even at high fre- 
quencies. 

Experimental confirmation of the ballistic w"^ regime 
is presented in Fig. [5l where the athermal component 
of the power spectral density (PSD) of both loosely and 
tightly crosslinked networks are given, as measured us- 
ing the same microrheology protocol as detailed in fl6{. 
Also shown are absolute shear moduli \G\ of both net- 
works over the same frequency range, confirming that the 
tightly crosslinked network is still in its elastic plateau 
up to frequencies of ~ lOOHz. Thus the downturn in 
the PSD is not due to the network response. It is how- 
ever consistent with the ~ w"^ ballistic regime predicted 
above, and we conclude this is its origin. Note that while 
there is evidence of a ~ regime at low frequencies for 
the sparsely crosslinked network, the data here is also 
consistent with a weak frequency variation, possibly due 
to the very slow, irreversible changes in network structure 
in response to strong force generations. 



III. ATHERMAL FLUCTUATION 
CORRELATIONS 

In this section, correlations in the fluctuations of the 
displacement field generated by a population of many 
firers are derived, for two classes of material. The starting 
point for the calculations is the displacement field around 
a single active firer, which will here be denoted by the 
material-dependent propagator P(x, t;n) that gives the 
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FIG. 5: (Color online) The upper panel shows the ather- 
mal component of the displacement power spectral density 
for loosely (filled squares) and a tightly (open triangles) 
crosslinked networks obtained from microrheology [l^ . The 
straight lines have the given slopes. The lower panel shows 
the absolute shear moduli |G(aj)| for the same networks, 
demonstrating plateau response and confirming the observed 
frequency-dependence in the power spectrum is not due to 
network viscoelasticity. 



displacement at a point x and time t due to the activity 
of a firer at x = and an earlier time t = 0. The total 
displacement u at a point r and time t can then be found 
by linear superposition of all active firers, 

u(r,t)=/dx/d./dnP(r-x,^-.;n)c(x,.;n) (4) 

where c(x, s; n) describes the number of firers per unit 
volume of orientation n at time s. We assume activity 
is statistically uniform and uncorrelated in both space 
and time, with a mean firing rate per unit volume (of 
any orientation) denoted c. Note that in an actin-myosin 
system, c refers to myosin mini-filaments and not simply 
the concentration of myosin. 



A. Mechanics of materials 

To proceed it is necessary to specify P(x, t; n) by de- 
scribing the mechanical response of the material. Two 
materials are considered here, both represented on the 
continuum level under the assumption that the quanti- 
ties of interest vary on lengths larger than the coarse- 
graining length of the network. First we consider the 
two-fluid model that has already been used to calcu- 
late the viscoelastic response of filament networks in so- 
lution [13, 113, HI]. In brief, the system is modeled as 
a coupled two-component system, consisting of a fixed 
isotropic elastic network interacting via local friction with 
an incompressible fluid. The network displacement u and 
solvent velocity v obey 



= MV^u + (^ + A)v(v-u) + r(v-c)tu) + r 
= y/V^v - VP - r(v - atu) 



(5) 



with /i, A the elastic Lame coefficients and rj the solvent 
viscosity. The pressure gradient VP imposes V • v = 
and the external force f describes all active force gen- 
erators. The frictional coupling coefRcient F is estimated 
to be ^ Tj/S,'^ with ^ the network mesh size jjlj. The 
solution is found by transforming to Fourier space, i.e. 



c(q,w;n)= J J dt e*("*-i-^)c(x, n) 



(6) 



and projecting the transformed u and v into components 
perpendicular and parallel to q, solvent incompressibility 
requiring that the latter vanishes for v. The transformed 
propagator P(q, w) can be written in terms of the pro- 
jected components of a single firer f^'^'(q, oj) as 

[P^*-q]q , P^^-[P^*-q]q 



where B — 2fi + X and G = fi — iurj is the total shear 
modulus. The force P^*(q, w) is the transform of either 
the oriented dipole ([1]), or the isotropic compression 



f-'(q,w) 



iK{uj)q[ci ■ n]n 
■|iK(a;)q 



oriented 
isotropic 



(8) 



The expression ^ is valid for F S> rjq^ 



9 < r 



which is anyway necessary to justify the use of the contin- 
uum equations [38]. The Poisson ratio i/ = X/[2{fj, + A)] 
can be estimated by assuming an affine network deforma- 
tion, giving the frequency-independent value = 1/4, as 
explained in Appendix |B] 

The second material considered here is a classical elas- 
tic body. Formally this corresponds to the F = limit 
in ([5|) , but allows the calculations to be carried out in real 
(rather than Fourier) space, which makes a small differ- 
ence to the calculation for single point fluctuations as ex- 
plained below. These calculations start from the known 
result for the displacement field due to a force monopole, 
from which the dipole field is found and the displace- 
ment correlations calculated by integrating the products 
of displacements. Details are given in Appendix [Cl 



B. Single— point athermal fluctuations 

Fluctuations in the network strain field can either be 
measured at the same point, or at points separated by a 
finite distance, corresponding to 1 and 2-point microrhe- 
ology, respectively. While fluctuations in active mate- 
rials have both a thermal and an athermal origin, they 
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may often be treated as independent and extracted sep- 
arately as described in |10|] . We calculate here the ather- 
mal contributions to these quantities deriving solely from 
firers, beginning with the 1-particle correlation function 
Cij{t) — {ui{O,0)uj{O,t)), where the angled brackets 
denote averaging over space, time and dipolc orienta- 
tions fi (isotropic compression centers are considered in 
Sec. IIII D|) . By making use of (jU and assuming uncorre- 
lated, uniform firers, (t) can be rewritten in frequency 
space as 




(a) 



(b) 



Cy(a;) = 



24^3 



dq|P(q,c^)|- 



(9) 



where (•••)& denotes averaging over orientation n only. 

For the propagators ([7]) and there is a short- 
wavelength divergence in Q at high q resulting from the 



slow decay of P, P 



The physical origin of this 



singularity is the diverging magnitude of displacements 
induced at positions arbitrarily close to an active firer. 
Recall however that the expression for P in ([7]) corre- 
sponds to point dipoles; for real finite-size dipoles, P will 
become strongly attenuated at length scales correspond- 
ing to the dipole size and this divergence will be avoided. 
In microrheology, a further small distance cut-off appears 
from the finite radius a of the probe particle, which ex- 
cludes dipoles from a range a of the origin. We focus on 
this form of cut-off and introduce a maximum magnitude 
for the wave vector, q"^^^ = n/2a and perform the 
integration in (j9l) using P defined by the oriented dipole 
in (O, giving 



607ra 



3|Gp 



= 1 - 



3? < Bz\/—i arctan 



(10) 



where the length £ — y^B/cuT is discussed below, 3?() de- 
notes the real part and the same root of should be 
chosen for both occurrences. Note that = \k{uj)\'^ is 
the power spectrum for a single firer, and there is implicit 
frequency dependence in both B{u}) and G(w). A repre- 
sentative example of (w) for an actin-myosin system 
is given in Fig. [TJ The corresponding result from the 
real-space elasticity calculations is 



Cf. 



CK 



(1 - i/)2 ISOtt^- 



■-Si 



(11) 



Note that this differs slightly from the low-frequency 

limit of m, = 0) = '^^j^^tMjJ^^S,,, by the 

constant coefficient in the numerator polynomial, which 
is attributable to the different manner of handling the 
short wavelength cut-off in real and Fourier space |31| . 



FIG. 6: Schematic representation of the effects of the cou- 
pling between the network and the solvent on the 1-particle 
response function Cijit) — {ui{O,0)uj{O,t)) . The crossover 
between the coupled response at large lengths (light grey) and 
the decoupled response at short lengths (dark grey) is given 
by the frequency-dependent penetration length £ = ^/B/uT, 
and the dominant response determined by the ratio of £ to 
the probe radius a. (a) For low frequencies, £/a ^ 1 and dj 
is dominated by decoupled network and solvent responses, 
(b) For high frequencies, £ <^ a and dj is now dominated by 
the coupled network-solvent response. 



Equation ([TU| is composed of shear /transverse modes 
oc |G|^2, and compressional/longitudinal modes oc |-B|~^. 
The contribution from transverse modes is independent 
of the network-solvent coupling, which is entirely con- 
tained in the longitudinal contribution. The magnitude 
of this term depends on the ratio of the cut-off rad ius a to 
the frequency-dependent penetration length £ = ^ B /ujT, 
which determines the range of the longitudinal network 
modes ^38]. Two different coupling regimes occur in the 
limits a <^ £ and a ^ as schematically shown in Fig. |6l 
For low frequencies such that a ^ £, the network is fully 
decoupled from the solvent and ip[z) 1. In this limit 
the longitudinal response is dominated by the elastic net- 
work. In the opposite limit a 3> network and solvent 
are fully coupled and now ?/)(z) w \B\'^tt^ / {%Qu:'^T'^a^). 



This term therefore decays as 



at high frequencies, 



so the longitudinal response vanishes and only the shear 
contributions remain. 



C. Two— point athermal correlations 

Spatial correlations of athermal fluctuations are 
encoded in the 2-particle function D^j (R, t) = 
(Mi(0, 0)m-, (R, t)), which quantifies correlated motion 
over a vector distance R. Using @ and again assuming 
spatio-temporally uniform activity, this can be rewritten 
in reciprocal space as 



(12) 



This can be evaluated following the procedure out- 
lined in Appendix |D] In terms of the pro- 
jections parallel and perpendicular to the line of 
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centers R, D^j{IV) = D^R)RiRj 
with R — |R|, we find 



D«{R,uj) 
D-^{R,Lu) 



D^{R){5,j - R,Rj) 



CK^ 




GOttR 1 








1207ri? 1 





3 
ISP 



^^(i?/^)^ (13) 



where = g?{2iB[l- (l + zV^)e-^^)]/z2}/s)fi(s) 

and ^Hz) = 3fJ{Be-^^^}/5R(B)-^^-L(z), and the root is 
chosen such that 5R(z-\/— i) > 0. See Fig.[l]for an example 
relevant for actin-myosin systems. For comparison, the 
real-space calculations give 



Df{R) 
Df{R) 



1 



607r/u2 R 



7 



1 



(1 



4807r/x2 R 



(14) 



which is exactly the zero-frequency limit of (|13p . 

Dij admits a similar scaling form as for Cy , but with 
the inter-particle spacing R instead of the probe radius a, 
Dij ~ cK-^/fM^R. As for dj, the coupling factors V'^ 
exhibit different high and low frequency behaviors, but 
now the crossover frequency is determined by the ratio 
R/£ instead of a/£, as schematically shown in Fig.[71 The 
decoupled limit arises at low frequencies such that £ ^ R^ 
when ip'^{z) — > 1 and ?/'" ~^ 0. Conversely in the fully 
coupled limit £ <C -R, both and -0" vanish. Note that 
0", ijj'^ are closely related to the response functions x", 
X"^ derived in [lllli. 

Despite the dependence of Dij on R through the cou- 
pling functions ■0''"'"(^/Oi only the overall scaling 1/R 
has been observed experimentally [1, [l^. The rea- 
son is simply that, for typical actin-myosin parame- 
ters, this additional source of variation is weak over the 
ranges of R typically probed, as confirmed in the inset 
of Fig. [1] The overall scaling R~^ is a direct conse- 
quence of the dipolar nature of the active agents. Di- 
mensional analysis suggests higher-order force multipoles 
will obey higher powers of R, i.e. for quadrupoles of the 
form /(r) ^ KQV^d{r) with kq the quadrupole moment, 
Dij (X cKq/fi'^R^ is expected (see Appendix [P]) . There- 
fore the experimental Dij ^ 1/R provides indepen- 
dent confirmation that activated myosin mini-filaments 
do indeed act as force dipoles to leading order. 

In the high frequency regime, or for an incompressible 
network with v = 1/2 for which B diverges, the coupling 
terms in (IT3t vanish and I?" ~ 2D-^. The same ratio be- 
tween parallel and perpendicular correlations also holds 
for thermal fluctuations in incompressible, passive me- 
dia 12|. Although one might hypothesize the existence 
of universality here, contrary evidence can be obtained 
by repeating the calculations for the force quadrupoles 
mentioned above. As described in Appendix |D1 in this 
case we expect Z?" = —2D-^, so the ratio between 




(a) 



(b) 



FIG. 7: Same as Fig. [6] but for the 2-particle correlation 
function _Dij(R,f) — {ui{O,0)uj{'R.,t)). The crossover be- 
tween the decoupled (a) and coupled (b) regimes is deter- 
mined by the inter-particle spacing R relative to the penetra- 
tion length £. Also shown in (b) are the regions responsible 
for negative and positive correlations to _D"''^°(7?,a;) when 
the firers are isotropic compression centers, as discussed in 

Sec. Imp] 



and D^ in athermal fluctuation is not solely determined 
by a material's mechanical properties. In particular, by 
mixing a population of dipoles and quadrupoles of suit- 
able moments it is hypothetically possible to generate a 
continuous range of ratios. 



D. Isotropic compression centers 



ISO 



For the isotropic compression centers given by ([2]) , Cj 
and Df° can be found by following the same procedure as 
above, with the propagator P"'°(q, w) in ([7]) now defined 
in terms of the isotropic force dipole from Repeating 
the calculations gives 



Df°{uj) 
DTicu) 



CK 



,5,,^(2a/TT£) 



1087ra|B|2 

,^Hr/£) 



CK 



36ttR\B\ 

CK^ 

72-kR\B\^ 



^^{R/£) 



(15) 



(16) 



where the coupling functions tp are the same as be- 
fore. Note that the shear terms have now vanished, 
and thus the variation in frequency due to coupling is 
far broader than before, suggesting a system of isotropic 
dipoles would present a stronger test for the prediction 
of the two-fluid model than the anisotropic case. Both 
Cij and Dij vanish in the strongly-coupled limit. Since 
■(/''^(z) — > as 2 — > 0, the parallel correlations w) 
also vanish at low frequencies, which is most easily ex- 
plained by reference to the real-space calculations below. 

Calculations for an elastic medium in real space are 
carried out following the same procedure as before, but 
for the displacement held due to an isotropic dipole. De- 
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tails are given in Appendix [C] The resulting expressions 
are 



D 
D 



isOjCl 

iso.cl 

II 

iso.cl 



1 - 2v 
l-v 



1 - 2v 
l-v 



CK 



4327r^2 Q 



1 



2887r^2 ji 



(17) 



(18) 



The results in all cases (including C-|°) match the low- 
frequency limits of the expressions |l5p and ([TCI) . As 
in ([T6)) . D^^° vanishes in the low- frequency, elastic limit. 
This is seen to arise, in real space calculation, from an 
exact cancellation of firers in between the two probe par- 
ticles, which generate negative correlations, and those far 
away, which generate positive correlations; see Fig. [71[b). 
This cancellation requires a displacement field that de- 
cays as the inverse-square with distance and would not 
occur for e.g. quadrupoles, which instead decay as 1/r'^. 



IV. DISCUSSION 

The expressions for athcrmal fiuctuations, Cij (uj) 
and DijCRjLu) derived in this study serve a number of 
uses. As explained above, the scaling of I?" and D-^ 
with distance R can be used to infer the leading-order 
force multipolc representative of the firers. In an active 
gel, myosin minifilaments might work either as oriented 
dipole or isotropic compression depending on their ag- 
gregation and geometries. By comparing the measured 
parallel and perpendicular correlations, it is also possible 
to determine the symmetry of the active firers (oriented 
or isotropic) . For materials obeying the 2-fluid model, we 
can further predict the extent of deviation from the in- 
compressible limits of Cij and Dij in terms of microscopic 
parameters, as expressed through the coupling functions 
-0, t/jII and tp'^ . In actin-myosin systems, this effect is 
found to be weak, accounting for a ~ 10 — 30% variation 
over the usual experimental frequency window, compara- 
ble to experimental precision. We can therefore extract 
the product c\k{uj)\'^ /\G{uj)\'^ , to good approximation, by 
fitting microrheology data to these expressions. Further 
measuring G{uj) independently with active microrheol- 
ogy, the quantity c|k(w)P can be determined. 

Active media are generally expected to have viscoelas- 
tic properties different from their equivalent passive sys- 
tems. For filament networks with active crosslinks {e.g. 
myosin mini-filaments), the forces generated by micro- 
scopic firers place the filaments in a state of non-zero 
mean tension, in contrast to a passive network which is on 
average unstrained under thermal fluctuations. Semifiex- 
ible polymers are known to stiffen under tension, thus an 
activated filament network is stiffer than a passive one in 
vitro [lO]. The local action of force firers c|«;(a;)p, which 
can be quantified using the expressions given in this 



study, is the microscopic origin of the observed macro- 
scopic active stiffening. Here we further provide two ten- 
tative scaling arguments that lead to expressions for the 
mean filament tension a in terms of c|K(a;)p, which can 
be used to predict the macroscopic pre-stressed network 
response [lO|. One argument is microscopic in nature, 
the second at the continuum level, and we show at the 
end are they in broad agreement with each other. 

For the microscopic argument, we first need an esti- 
mate of the number of active myosin crosslinks per actin 
filament. Taking the volume associated with a filament 
of length £ in a network of mesh size ^ to be ~ ££,'^, 
then the typical number of myosin crosslinks per filament 
should scale as ~ c^^^. Substituting typical values for in 
vitro actin-myosin systems (^ ~ lOOnm, £ ~ 1 — 10/im 
and c ~ l^m~^ (67]]) gives < 0(1) active crosslinks per 
filament. We therefore ignore the possibility of multi- 
ple active forces on the same filament. The longitudinal 
component of a force applied to a filament will propa- 
gate a distance ^ A before being transferred into bending 
modes of filaments in the surrounding network, where A 
is related to the mean distance between crosslinks ic (as 
measured along the filament contour) by A ~ £c{ic/dA)^ 
with z a small exponent and cIa the characteristic fila- 
ment diameter [4^, . Therefore the mean tension a 
in the network is estimated to be the number of active 
crosslinks per filament scaled by X/£, multiplied by the 
force per active crosslink ^ k/e. Assuming that the mo- 
tor generated forces act between crosslinks [Toj, e ^ £c: 
the mean tension in the active network is obtained as 



(19) 



The exponent z has been predicted to be 2/5 [4a], al- 
though the slightly different value z = 1/3 was empir- 
ically better in athermal, two-dimensional networks at 
intermediate crosslinking densities psj . 

For the continuum calculation, we consider a single 
oriented force dipole in an isotropic elastic body, which, 
when active, generates the strain field Uij derived in Ap- 
pendix [Cl The relative extension of a filament segment 
of length ic is estimated as Sl^jlc 

Hi 1J Til 2 Thj with n the 

filament orientation |68j , assuming Uij is uniform over its 
length. Averaging Uij over all filament and dipole ori- 
entations trivially gives a zero mean extension 5£c — 0. 
However, filaments can only support limited compression 
Sic < before relaxing their longitudinal stress via buck- 
ling. Therefore the integration should be restricted to 
purely elongating filaments 6£c > within a range X 
from the dipole. At this length X, the strain is just large 
enough to initiate buckling; beyond X, both compres- 
sion and elongation are allowed and the contribution to 
S£c averages to zero. Averaging and integrating Uij with 
this restriction, we obtain that the mean extension per 
filament is S£c/£c ~ (ck / n) \n{X / £c) . The logarithmic 
factor \n{X/£c) derives from the inverse-cube decay of 



the strain field with distance r, Ui 



, necessitatmg 
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both long and short wavelength cut-offs. Since it is log- 
arithmic, it varies only slowly with X and £c and may 
be regarded as a constant in most cases. To calculate 
the mean tension a, we use the known force-extension 
relation a ~ iilk^^£^){Sic/Q with ip the filament's 
persistence length [39|, giving the final estimate 



e;kT CK 



ln(X/4) 



(20) 



Regarding the weakly-varying multipliers ln{X/£c) and 
{lc/dA.Y^^ as numerical prefactors, (fT9|) and ([20l) give 
us the network shear modulus /x ~ l^k^T / i^^'^ in agree- 
ment with prior studies [39| . Despite the very differ- 
ent approximations taken, these two estimates are thus 
broadly consistent. Therefore by determining c|k(cli)P 
and G{lo) using microrheology experiments and the ex- 
pressions derived in this manuscript, it is now possible to 
completely address the mechanism of how the nonequi- 
librium activity leads to the modification of the network 
viscoelasticity via (|19p or pO|) . 



first term on the right hand side of (jAip corresponds 
to an initially growing k that does not start releasing 
by time t, the second to those that start releasing at 
some time < s < t and are still in this release phase 
at t, the third to initially releasing k, and the fourth 
and fifth terms to second and subsequent cycles after 
the first release phase has finished. The derivation of 
(|Aip is straightforward; for instance, the second term is 
Pa /j^ ds ie^^/'^(5(K — [kq + kqs — (t — s)kr]), where the 
exponential gives the probability for the first release 
event at s, the delta function maps to the corresponding 
value of K at t, and ti = max(0, [tkj^ — ko]/[^g + ^fl]) 
ensures k > 0. Performing the integral gives the second 
term in (IA1[) . The fourth and fifth terms are derived in 
a similar manner. We note that Q{K,t;KQ) is properly 
normalized, dKQ{K,t; kq) = 1 Vt, kq. 

As i ^ oo, Q{K,t;Ko) -> Qoo(k) = T^^e-''/^''^'^^ .^rhich 
is just an exponential distribution with mean (k) = tkq. 
The correlation function C{t) — (^(io + ^)'^(^o))to ~ (^)^ 
in steady state is found by integrating over (jAip suitably 
weighted by Qoo(ko), 
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C{t) = / dKo / dKKoKQoo(Ko)Q(K,i; Ko) - (k)^ 
Jo Jo 



1 



1 



1 - a 



1 — a 



(A2) 



Appendix A: Finite release— time model 

Here we derive the power spectrum for the 

stochastic grow-and-release model given in Fig. [31 In 
steady state, the dipole moment /t(t) will be growing 
with probability = kr/ (kq + /^r) and releasing with 
probability pR = 1 — PG = kG/{ko + kr). To find 
Q{K,t;uo), the probability density function of k at time 
t given AC = at t = 0, we therefore integrate over all 
possible histories and weight according to pq or pR. The 
result is 



= pGe~*^^S{K- [tkc + Ko]) 



using a — kq/kr < 1 for the ratio of growth to release 
rates. This exhibits positive correlations on the charac- 
teristic time scale for growth r, and negative correlations 
on the release time scale ar. In the rapid release limit 
a — we recover the expected correlation for a Poisson 
process, C{t) = (K)2e-*/^. Note that although we as- 
sume a < 1 for this discussion and in the main text, the 
result is also valid for a > 1 and indeed the steady state 
is invariant under {a, kg, kr, t) ^ i<^^^, kr, kg, cxt). 

It is now straightforward to Fourier transform (jA2p 
to get the power spectrum for a single firer — 
J^ao di e*"*C(|i|), which is given in ([S]). For completeness 
we also give here the limiting case kq — > kr when the cy- 
cle is symmetric. We find that C{t) ^ t/r)e~*''^ 
and |K(a;)p = 4:t{k)'^ / , so the pseudo-diflFusive 
regime is now entirely absent. 



-1- pGd{K — [ko — t^is])(^{no + tko — n)e 



Tkr 



+ pr(5(k - [ko - ^kr]) 

+ p^e{t - ko/kr)(5(k -kG[t- Ko/KR])e-*/"e«''/"'^'' 



KR KG 



1 -k/tkg -t/T 1 

- — e ' — Pgg — 



TKQ 



Tkji 



(Al) 



where T — r(KG + 'iR)/'^R is the mean total cycle 
time (growth plus release), 6 is the Heaviside step func- 
tion and K > is implicitly assumed throughout. The 



Appendix B: Poisson ratio for afflne network 
deformation 

Here we derive the Poisson ratio for a 3-dimensional 
filament network under the assumptions of affine quasi- 
static deformation and a homogeneous, isotropic net- 
work. Consider a single filament of length £ and ori- 
entation n with respect to some fixed axes. If the system 
is deformed to give a displacement field u{x,t), and as- 
suming this macroscopic field scales uniformly down to 
lengths <C £ {i.e. the deformation is affine), then the 
filament elongates by an amount S£ given by [6a | 
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£ — UijUiTlj 



(Bl) 



where Uij = (diUj + djUi)/2 is the linear strain tensor. 
The increase in energy in this filament's deformed state 
is then 5E = ^ki^Stf' = ^U'^{uijnihjY with k an effec- 
tive spring constant. Note that k may depend on ^, the 
frequency of deformation change, any pre-tensile stress 
in the filament etc. For the purposes of this calculation, 
however, we need to only assume k is independent of the 
filament orientation n. 

Now consider a homogeneous population of filaments 
with length distribution C{t), orientation distribution 
A/'(n) and number density p. To ensure isotropy, t and 
n are uncorrelated and A/'(n) — (47r)~^. The increase in 
energy per unit volume is then 



TABLE I: The displacement, strain and stress fields for an 
isotropic compression center of size e centered on the origin. 

g(5ij — fifj. 



Field 


r < e/2 


r > e/2 


Displacement uf°(r) 


-2v) r . 


k(1 -2u) 1 . 


37r/i(l ~ v) 


247r^(l - v)r'^ 


Strain «Jj°(r) 


k{1 -2v) 1 


n[l-2v) 1 ^ 




87r/i(l — v) '■^ 


Trace uf°{v) 


k(1~2v) 1 
7r^(l — v) 





Stress of°{v) 


2k{1 + u) 1 


41 "2i/) 1 ^ 


37r(l - I.) ' 


47r(l — v) 



dlC{l) j dn7V(n)ifc£2(u„n,rij)2 (B2) 
dhiuijhifijY (B3) 



<r(r;fi) 



47r 



where T, = J d£C{l)k£^. Assuming quasi-static re- 
sponse, this energy change equals the corresponding elas- 
tic energy Ai?°' = ^Xuf^+pufj with A, p the Lame coeffi- 
cients [68] . By choosing two independent strain fields Uij 
it is therefore possible to determine both A and p. For 
instance, an isotropic compression Uij — —jSij gives the 
equation 7^1] = |A7^ -I- Spj'^. For simple shear, where 
= 7/2 and Uij = otherwise, another equa- 
tion 7^I]/15 = Ai7^/2 is obtained. From these 2 equa- 
tions, we obtain X ^ p = 2S/15, and the Poisson ratio 
= A/[2(A + p)] = 1/4, independent of S and hence k. 
By coincidence, this is in good agreement with the value 
obtained for semi-dilute flexible polymers in a good sol- 
vent 



Appendix C: Real-space calculations 

The displacement field u(r) due to a point 
force f applied at the origin of an infinite, ho- 
mogeneous, isotropic elastic continuum is known, 
u^{r) = [167r^(l - vy]"^ {(3 - 4t/)(5y + fifj} fj, with p 
the shear modulus, v the Poisson ratio, r = |r| and 
r = r/r The displacement field extending from the 
dipole with orientation n defined in ^ is then obtained 
by taking the gradient of the monopole solution. 



(r;n) 



167r/i(l — ;/) r 
-(l-3[f.np)f,} 



(CI) 



with K ~ ef the dipole moment. The linear strain Uij = 
{diUj + djUi)/2 and stress fields can be found by the usual 
relations. For instance, the strain is 



167r/i(l — 1/) 
|2(1 - 2i')fiihj - 3(5[f • n]2 - l)f,fj 

+ (3[f • n]2 - l)Sij + 6i/[f • h]{nifj + (C2) 

This equation has a trace m,™ = — k(1 — 2z^)[87r/i(l — 
iy)r^]^^{l — 3[f • n]^) that vanishes in the incompressible 
limit J/ = 1/2 as expected. 

It is possible to derive an exact solution for finite dipole 
with size e by exploiting the spherical symmetry that 
arises after averaging over n. The resultant displacement 
field is found by integrating the monopole solution over 
a sphere of inward-pointing monopole forces — 2/f . The 
sphere has radius s/2 centered on the origin, and the fac- 
tor 2 arises since each dipole consists of 2 monopoles. By 
symmetry only the response at, say, a point (0, 0, z) needs 
to be considered. The integration is simplified by switch- 
ing to spherical polar coordinates and making the change 
of variables u — yjr^ \ e^/4 — re cos 9 for the polar angle 
6; the cases z < e/2 and z > e/2 correspond to different 
integral limits and must be considered separately. The 
results are presented in Table HI 

The firers are assumed to be independently and uni- 
formly distributed throughout space, so the firer concen- 
tration field c(x) obeys (c(x)c(y)) — c^(x— y) and hence 



= Qd^J dyc(x)c(y)7.,(~x)7., (R-y)^ (C3) 

= dxwi(-x)Mj (R - x) (C4) 

The single-point correlations correspond to R = 0. The 
required expression for a single isotropic or anisotropic 
field Ui is then inserted into the integrand. For the 2- 
point correlations it is convenient to set R = (0,0, i?) 
for the integration, and rotate to general R at the end. 
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The integration is facilitated by switching to spheri- 
cal polar coordinates by converting variables as u = 
y/r^ + E? ~ 2rR cos 9. We give here an intermediate ex- 
pression for the correlations due to a nematic field of 
anisotropic dipoles with the same orientation n, 



7-jani 



CK 



;(l-[n.R]2)x 



|(1 - [n • R]^)(5jj + 4[n • R](2 - 3iy)[n,Rj + R^fij] 
+4(1 - 2i/)(3 - ^v)h.,nj - (1 - 5[n • AY)RiRj^ 

(C5) 

The expression given in (I14[) is then found by averaging 
(C5) over n. For isotropic dipoles we assume e/r 
and hence only use the r > e/2 expressions in Table ID 

Appendix D: Dij for the 2— fluid model 

Dij (q, Lj) for the 2-fluid model is found by inserting the 
anisotropic dipole propagator P(q, from (O and (|8]) 
into (jl2p and then averaging over orientation n. In terms 
of projection operators parallel and perpendicular to q, 



= ck' 



q 



1 



15(72|G|2 



(Dl) 



Note this projection is in q-space, not real space. 
For an isotropic material £'y(R) can similarly 



be projected into components parallel and per- 
pendicular to the real-space unit vector R, i.e. 
Aj(R) = D\\{R)R,Rj + D^{R){S,, - R,Rj). 

To relate the q-space projections to the real 
space ones, note that RiRjDijiTL) — D^^{R) and 
((5y - R,Rj)D,j{IC) = 2D^{R). By writing Aj(R) as 
the inverse transform of (|D1I) in both of these expres- 
sions, it is readily seen that the following transform pairs 
hold between ^"'-^(i?) and -Ell'^(g), 



£;^(g)(l-[q.R]2)+i?ll(g)[q.R]2 
£;^(g)(l + [q.R]2) + i?ll(q)(l-[q.R]2) 

(D2) 



Z?"(i?) ^ 
2D^{R) ^ 



Performing the inverse transforms of the E^^'-^{q) in (jDl[) 
using (|D2[) gives the final expressions (|13p . The inte- 
grations are facilitated by switching to cylindrical po- 
lar coordinates {p, 9, z) and exploiting isotropy to set 
R = (0, 0, R) without loss of generality. 

It is possible to make some general remarks on 
the correlations of displacement fields driven by force 
quadrupoles, /(r) ^ kqW^5{v) with kq the quadrupole 
moment, without going through detailed calculations. 
After averaging over the quadrupole internal degrees of 
freedom, Dij is expected to take the same form as IjDip 
but with E^{q) cx CKg/jCp, the different power of q 
(with respect to dipoles) being due to the one-higher 
derivative of 5{v) in f. The longitudinal component 
= for an incompressible body. Repeating the in- 
verse transformations (jD2[) now gives oc 
ckq/\G\'^R^ , as discussed in the text. 
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